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ABSTRACT 

We use deep HST ACS/HRC observations of a field within M32 (Fl) and an M31 background field (F2) to 
determine the star formation history (SFH) of M32 from its resolved stellar population. We find that 2-5 Gyr 
old stars contribute ^ 40% ± 17% of M32's mass, while ^ 55% ± 21% of M32's mass comes from stars older 
than 5 Gyr. The mass-weighted mean age and metallicity of M32 at Fl are (Age) = 6.8 ± 1 .5 Gyr and ([M/H]) = 
-0.01 iO.OSdex. The SFH additionally indicates the presence of young (< 2 Gyr old), metal-poor ([M/H] ~ 
-0.7) stars, suggesting that blue straggler stars (BSS) contribute ^ 2% of the mass at Fl; the remaining ^ 
3% of the mass is in young metal-rich stars. Line-strength indices computed from the SFH imply a light- 
weighted mean age and metallicity of 4.9 Gyr and [M/H] = -0. 12 dex, and single-stellar-population-equivalent 
parameters of 2.9 ±0.2 Gyr and [M/H] = 0.02 ±0.01 dex at Fl (^ 2.1 re). This contradicts spectroscopic 
studies that show a steep age gradient from M32's center to 1 r,,. 

The inferred SFH of the M31 background field F2 reveals that the majority of its stars are old, with ^ 95% 
of its mass already acquired 5-14 Gyr ago. It is composed of two dominant populations; ^ 30% ±7.5% of 
its mass is in a 5-8 Gyr old population, and ^ 65% ± 9% of the mass is in a 8-14 Gyr old population. The 
mass-weighted mean age and metallicity of F2 are (Age) = 9.2± 1.2Gyr and ([M/H]) = -0.10 ± O.lOdex, 
respectively. Our results suggest that the inner disk and spheroid populations of M3 1 are indistinguishable 
from those of the outer disk and spheroid. Assuming that the mean age of M3rs disk at F2 (^ 1 disk scale 
length) to be ~ 5-9 Gyr, our results agree with an inside-out disk formation scenario for M3rs disk. 
Subject headings: Local Group — galaxies: individual: M32, M31 — galaxies: elliptical and lenticular, cD — 
galaxies: star formation history, stellar content 



1 . A BRIEF fflSTORY OF STAR FORMATION IN M32 

M32 (NGC 221) is a compact, low-luminosity elliptical 
galaxy, satellite of our neighbor M31. Due to its proximity. 
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we can study M32 with great detail not only from its inte- 
grated light but also from its individual, resolved stars in a 
way that is impossible for most of the elliptical galaxies, given 
their greater distances and high densities. Thus, M32 is a very 
important galaxy to understand the formation and evolution 
of low-luminosity spheroidal star systems. However, M32's 
SFH, and therefore its origin, is still controversial. The dif- 
ferent scenarios proposed to explain its origins extend from 
a true elliptical galaxy at the lower extreme of the mass se- 



quence (e.g., Faber|197"3 Nieto & Prugniel|1987[ Kormendy 
et al. 2009 [ l to a threshed spiral galaxy (e.g., Bekki et al. 2001 
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Fig. 1 .— Location of our two HST ACS/HRC pointings: M32 (Fl) field 
and M31 background (F2) field, indicated as small black boxes. Each field 
covers a region of 26 X 29arcsec^ on the sky. The field Fl is located at 
1 10" from the nucleus of M32 and represents the best compromise between 
minimizing image crowding and contamination from M31. The F2 field is 
at the same isophotal level in M31 as Fl. At the distance of M32. each field 
occupies an area of 11752 pc^. Thirty-two exposures in each of the F435W 
(B) and F555W (V) filters were taken for each field. North is up and East is 
to the left. 



[Chilingarian et al.'2009'). 

The only way to accurately determine the age, and thus the 
SFH, of a galaxy is by directly observing its oldest main- 
sequence tumoff (MSTO). With this goal in mind, we were 
awarded 64 orbits with HST ACS/HRC to observe two fields 
near M32, Fl and F2 (Fig. [TJ, in order to detect the oldest 
MSTOs of this galaxy. 

1.1. The deepest HST CMD ofM32 
In |Monachesi et al. (2011 hereafter Paper I) we intro- 



duced our observations and presented the deepest HST color- 
magnitude diagram (CMD) of M32 yet obtained, reaching 
more than 2 mag fainter than the RC and fully resolving 
the RGB and the AGB. Paper I significantly improved our 
knowledge on the stellar populations of M32. We have found 
that M32 is dominated by intermediate-age (2-8 Gyr old) 

and old (8-10 Gyr old), metal-rich ([Fe/H] 0.2) stars 

and it contains some ancient (> 10 Gyr), metal -poor stars 
([Fe/H] ^ -1.6) as well as possible young populations (0.5 
- 2 Gyr old stars). 

These conclusions were provided by our qualitative analy- 
sis of the CMD of M32, which shows a red clump (RC), a red 
giant branch (RGB), a RGB bump (RGBb), an AGB bump 
(AGBb), and a blue plume (BP). Figure 12 of Paper I, repro- 
duced here as Figure [2] shows a Hess representation of the 
CMD of M32 decontaminated of M3 1 stars, where the differ- 
ent evolutionary features are highlighted. We summarize here 
the main findings and conclusions of Paper 1. 

• The core-helium burning stars are concentrated in a RC 
and its mean color and magnitude suggest a mean age 
of 8-10 Gyr for a metallicity of [M/H] - -0.2 in M32. 

• The first detection of the RGB bump and the AGB 
bump in M32 permits a constraint on the mean age and 
metallicity of the population. This gives a mean metal- 
licity of M32 higher than [M/H] ^ -0.4 dex and a mean 
age between 5 and 10 Gyr. 

• The metallicity distribution of M32 inferred from the 
CMD has a peak at [M/H] ~ -0.2dex. Overall, the 



metallicity distribution function implies that there are 
more metal-rich stars than metal-poor ones. Metal- 
poor stars with [M/H] < -1.2 contribute very little, at 
most 6% of the total V-light or 4.5% of the total mass, 
to M32 in Fl, implying that the enrichment process 
largely avoided the metal-poor stage. 

• Bright AGB stars at F555W < 24, i.e. above the TRGB, 
confirm the presence of an intermediate-age population 
in M32 (ages of 1-7 Gyr). 

• The observed blue plume is genuine, not an artifact of 
crowding, and contains stars as young as ^ 0.5 Gyr. 
The detected blue loop, with stars having masses of 
~ 2-3 M0 and ages between ^ 0.3 and ~ 1 Gyr, and 
the possible presence of a bright SGB are different man- 
ifestations of the presence of a young population. How- 
ever, in Paper I we suggest that it is likely that this 
young population belongs to the disk of M31 rather 
than to M32. The fainter portion of the blue plume 
(F5551V > 26) does belong to M32 and indicates the 
presence of stars with ages 1-2 Gyr and/or the first di- 
rect evidence of blue straggler stars (BSS) in M32. 

• The oldest MSTOs were out of reach, given the severe 
crowding in Fl, and there is no significant BHB ob- 
served in Fl, so an ancient, metal-poor population can- 
not be seen directly in our CMD. We have, however, a 
hint of the presence of such a population from a l-a 
detection of RR Lyrae stars found in Fl and associated 
with M32 using our data ( [Fiorentino et al.|20lO] l. 

• In general the CMDs of both fields Fl and F2 show an 
unexpectedly similar morphology. By subtracting the 
normalized Fl CMD from the F2 one (see Figure 21 
in Paper I), one can detect subtle differences. M31 has 
a younger and more metal-poor population than M32, 
and M32 has a more conspicuous intermediate-age pop- 
ulation (Fig.|4]i. 

• The CMD of our M31 background field F2 exhibits a 
wide RGB, indicative of a metallicity spread with its 
peak at [M/H] ^ -0.4 dex. The presence of a blue 
plume indicates the presence of stars as young as 0.3 
Gyr. Bright AGB stars in F2 reveal the presence of an 
intermediate-age population in M3 1 . 

1.2. Completing the picture ofM32's SFH 

The analysis presented in Paper I provided initial con- 
straints on the ages and metallicities of the stellar populations 
of M32 at Fl and M31 at F2. That work was based on tradi- 
tional methods of isochrone analysis, and was heavily based 
on the brighter evolved portions of the CMDs, such as the RC, 
RGB, and bump (RGBb and AGBb) features. 

The approach in this paper is independent. Here we use 
a sophisticated method of CMD analysis and decomposition 
that digs into the fainter, severely-crowded portions of the 
CMDs near the MSTO and sub-giant branches. We recover 
information from the brighter MSTOs present in the CMDs, 
thus providing quantitative information about the younger 
populations of M32. In this paper, we derive the detailed 
young and intermediate-age SFH of M32 at ^ 2' from its cen- 
ter and of M31 at our background field's location, which was 
not possible from the analysis in Paper I. We find that our field 
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Fig. 2. — Error-based Hess diagram for M32, corrected for contamination from M31 background stars. The boxes indicate various features that represent 
different stellar populations. MS: Main Sequence; BP: Blue Plume; SGB: Subgiant branch; BHB: Blue Horizontal Branch; BL: Blue Loop; RC: Red Clump; 
RGBb: Red Giant Branch bump; R-RGB: Red-Red Giant Branch; B-RGB: Blue-Red Giant Branch; TRGB: Tip of the Red Giant Branch; AGB: Asymptotic 
Giant Branch; and AGBb: Asymptotic Giant Branch bump. The dotted-dashed line indicates the 50% completeness level of our data. Magnitudes are calibrated 
onto the VEGAmag system. (This is Fig. 12 from Paper I; we refer the reader to that paper for more details.) 



in M32 has a substantial population of 2-5 Gyr old stars con- 
tributing to ^ 40% ± 17% of its mass, an unexpectedly large 
population of young stars at such a large distance from the 
center of an elliptical galaxy. 

The paper is organized as follows. In Section |2] we briefly 
describe our observations and photometry. Section [3] de- 
scribes the method used to derive the SFH. We present the 
results of the SFH analysis obtained for Fl, F2 and M32 in 
Section |4] In Section |5] we provide a detailed and complete 
SFH of M32 and discuss its implications on M32's origins, 
synthesizing a complete picture based on both the present and 
Paper I analyses. In Section|6]we discuss the SFH of the in- 
ner regions of M3 1 . Finally, we summarize our results and 
present our conclusions in Section]?] 

2. OBSERVATIONS AND PHOTOMETRY 

The field selection and observational strategy as well as the 
image reduction are described in Paper I and we refer the 
reader to that paper for details. Briefly, HST ACS/HRC im- 
ages of two fields near M32 were observed during Cycle 14 
(Program GO-10572, PI: Lauer). The M32 HRC field (Fl) 
was centered on a location 110" south (the anti-M31 direc- 
tion) of the M32 nucleus. The background field (F2) was lo- 
cated 327" from the M32 nucleus, roughly along its minor 



axis, at the same isophotal level in M31 as Fl. The field loca- 
tions are shown in Figure [T] 

Each field was observedror 16 orbits in each of the F4-35W 
(^ B) and F555W (^ V) filters. All of the images were 
combined in an iterative procedure designed to detect and re- 
pair cosmic -ray events, hot pixels, and other defects, with a 
Nyqu ist-sampled summed image as the final product (l Lauer| 
|1999| l. Color images of Fl and F2 are shown in the left and 
right panels of Figure [3] respectively, where the strong crowd- 
ing in these fields is clearly visible. There is however a differ- 
ence between the stellar density in Fl and F2: the crowding 
is more severe in Fl than in F2. This can also be seen from 
the bottom panels of Figure[3] where zoomed-in images of the 
top panels are shown. 

Stellar photometry was performed on deconvolved com- 
bined images. A detailed description of the deconvolution 
process is explained in Paper I. In short, deconvolution was 
performed on t he final ima ges using the Lucy-Richardson al- 
gorithm (Lu cy J974[ [Richardson 1972) and empirically con- 
structed PSFs, one for each image. Stars were identified in the 
deconvolved images and their fluxes were measured. Change 
transfer efficiency (CTE) and aperture corrections were ap- 
plied to the magnitudes, which transform the instrumental 
magnitudes into calibrated, apparent magnitudes. 
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Fig. 3 . — Combined color images of tlie 32 exposures in the Fl (top left panel) and F2 (top right panel) fields displayed with the same logarithmic stretch. Each 
image has a size of 2048x2048 pixels with a 0"0125 pixel scale. There is a clear difference in stellar density between the images, indicating that crowding is 
more severe in Fl than in F2. We also note a stellar density gradient in the Fl image, becoming higher when approaching the center of M32. The long black spot 
in the top center of each image is the occulting finger of the ACS/HRC coronagraph. The bottom panels are zoomed-in images of the centers of the top images 
for Fl (left) and F2 (right) fields, where we can better see individual stars and the different crowding levels are also evident. Note the blue stars in these images. 
Each zoomed-in image represents an area of ~ 15arcsec^ on the sky. 



Figure |4] shows the CMDs derived for Fl (left panel) and 
F2 (right panel) from the deconvolved photometry, calibrated 
onto the VEGAmag system. They contain 58143 and 27963 
stars, respectively, as indicated in Table [T| A qualitative anal- 
ysis of these CMDs allowed us to gain some insights into its 
stellar populations. This was discussed in detail in Paper I and 
we have summarized our conclusions above. 

Note the difference between the CMD of Fl and F2 at mag- 
nitudes between ^ 27 and 28 (cyan boxes in Fig- 
ure |4|l . The number of stars in this region, where the brighter 
MSTOs are located, is larger in Fl than F2. This suggests that 



there is a bigger contribution of intermediate-age stars in Fl 
than in F2. We can better appreciate this difference in Fig- 
ure 21 of Paper 1, where we showed a Hess subtraction of the 
normaUzed Fl CMD to the F2 CMD. 

2.1. Crowding tests 

We performed artificial star tests (ASTs) to assess the com- 
pleteness level and quantify the photometric errors of our data. 
This is a crucial step for the derivation of the SFH. The distri- 
bution of stars in the observed CMD is modified from the ac- 
tual distribution due to the observational effects, particularly 
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TABLE 1 
Deconvolved photometry 



Field 


Detections" 


nF*35Wb 
'^PSF 


nf5S5Wb 


ACf435w'- 




Fl 

F2 


58,143 
27,963 


5 
6 


5 

16 


-0.25 
-0.22 


-0.22 
-0.10 



Final number of stars detected and used to derive CMDs 
^ PSF radius in HRC original pixels 
^ Aperture correction 

at the fainter magnitudes where most of the information from 
the older star formation is encoded. The ASTs are used to 
simulate the observational effects in the synthetic CMDs that 
are then compared with the observed CMDs in the analysis 
described below. 

The procedure and results of the ASTs are presented in 
Paper I and we refer to that paper for further details; we 
give a brief description here in order to provide g uidance for 
later sections of this paper. We used lAC-STAR (Aparicio & 



[Oallart 200 4 ) to generate 5 x 10^ artificial stars with realis- 
tic colors and magnitudes covering not only the entire color 
and magnitude range of the observed stars but also ^ 2 mag- 
nitudes fainter We injected the artificial stars into the real 
images after transforming their magnitudes into instrumental 
ACS/HRC fluxes. The number of stars injected per experi- 
ment is 2000, to avoid increasing the already severe crowding 
of the real images. We performed 250 ASTs per field/filter 
combination for a total of 1000 ASTs. The images containing 
real and artificial stars are photometered exactly in the same 
way as the original images. A comparison of the known in- 
jected magnitudes and colors of the artificial stars to those 
obtained from their photometry allows us to quantify the pho- 
tometric errors. The completeness of our data at a given 
color and magnitude is calculated as the ratio of recovered- 
to-injected artificial stars on that color and magnitude bin. 

The results obtained from these ASTs indicate that the lim- 
iting magnitudes of the Fl and F2 CMDs are F555W 28 
and ^ 28.5, respectively, nearly independent of color. The 
CMD of F2 is therefore slightly deeper than that of Fl (cf. 
Figs. 8 and 9 of Paper I). The 50% completeness level as well 
as the photometric errors derived from the ASTs for Fl and 
F2 are indicated in Figure |4] 

3. THE lAC METHOD TO RESOLVE THE SFH 

To extract the detailed SFH of Fl and F2 we use the well- 
known met hod of fitting synthetic CMD s to the da ta (see e.g., 
Tosi et al.|[ l991; BertelH et al. 1992; Tolstoy & Saha 1996; 
Aparicio et al. 1997). There are currently many approaches 
to derive detailed SFH of galaxies (e.g., StarFISH: Harris 



Zaritsky 2 00T] MATCH: |polphin||2002l lAC- pop/f^innlAC: 
Aparicio & Hidalgo|2009 |Hidalgo e t al. '20111 as well as sev- 
eral d ifferent stellar libr anes (e.g., B aSTI: Pietrinfemi et al. 
SOOT Padova/Girardi: ' Girardi et al.||2000f jMarigo et al. 
^08 ) available to compute the required synthetic CMDs. We 
use the lAC-pop/MinnlAC method and adopt the Ba STI and 
Padova stel lar libraries. The lAC-pop code (Aparicio & HiJ 
dalgo 2009 1 uses a modified merit-function (Mighell 1999JI 



to compare the observed and synthetic star counts in different 
boxes (see bel ow) of the CMDs. A genetic algorithm ( Char^ 
bonneau 1995 1 is adopted to minimize x^- An important char- 



acteristic of the code is that it solves the SFH simultaneously 
for age and metallicity distributions. It thus provides the SFH 
of a stellar system as a linear combination of simple popula- 
tions, i.e. within small ranges of age and metallicity. We refer 



the rea der to Aparicio & Hidalgo ( 2009| l and Hidalgo et al. 
( |201 1| ) for more details. 

It is important to emphasize that, for the current analy- 
sis, we have mainly used information from the extended MS, 
MSTO and SGB regions of the CMDs, as we will see be- 
low. We have excluded the RC and most of the RGB regions, 
which were the main features analyzed in Paper I and from 
which we obtained estimates on the age and metallicity of 
M32. This is because the physics governing these phases are 
more uncertain than those on the MS and SG B, and differ- 
ences between stellar libraries are more severe (|Gallart et al.| 



2005j l. For instance, the morphology and number of stars oc- 



cupymg the HB/RC evolutionary phases depend on unknown 
issues, like mass loss on the RGB or He-core mass. Small 
differences in the adopted physics can significantly alter the 
number of stars and morphology of these CMD regions. The 
CMD regions that we probe in this paper allow us to obtain de- 
tailed information about the young and intermediate-age pop- 
ulations of M32, something that was not possible in Paper I, 
but conversely, we cannot make a quantitative analysis of the 
older populations and so we must rely on the qualitative re- 
sults of Paper I. 

3.1. Steps carried out to obtain the SFH 

1) Synthetic C MD. We first generate a s ynthetic CMD us- 
ing lAC-STAR (j Aparicio & Gallart|[2004| . The bolometric 
corrections applied to both libraries are those of |Origlia &| 
Leitherer (2000) which transform the theoretical tracks into 
the ACS/HRC photometric system. We assume a constant star 
formation rate (SFR) from to 14 Gyr, and metallicities from 
Z = 0.0001 ([M/H] = -2.3) to Z = 0.04 ([M/H] = 0.3) uni- 
formly distributed at all ages. Note that there is no assumed 
age-metallicity relation as input, and the selected age and 
[M/H] ranges are broader than those expected for the solu- 
tion. This allows the code to find the SFH solution with mini- 
mum cons t raints and ensures no lost information. We adopted 
a |Kroupa| ( |2002| initial mass function (IMF)^ from 0. 1 to 
1OOM0. The IMF has a slope of 1.3 for stars with masses 
lower than 0.5 Mq and 2.3 for stars with higher masses. We 
assume a 35% binary fraction with a relative mass ratio ran- 
domly distributed between 0.5 and 1 (the impact of differ- 
ent binary fractions on the solution is discussed in the Ap- 
pendix). The synthesized CMD, shown in the left panel of 
Figure pi contains 5 x 10^ stars and its faintest magnitude is 
^ 2 magnitudes fainter than the 50% completeness level of 
our data. Observational effects (incompleteness and photo- 
metric errors) are simulated using i nformation obta ined from 
the ASTs described above (see Hida lgo et al.|201 1 and refer- 
ences therein for a detailed description of this procedure). The 
right panel of Figure |5] shows the synthetic CMD after obser- 
vational e ffects are simulated. We call it a "model CMD" 
following | Aparicio et al.| ( [T997 )'s notation. The model CMD 
is the one to be compared with the observed CMD for the 
derivation of the SFH of our fields. 

2) Parametrization of the CMDs. This is the main input 



of the lAC-p op code and was performed using MinnlAC ( Hi- 
|dalgo et aLpOl 1 ), a set of routines specially designed for this 
purpose. We first define the "simple populations", the age 

' Given that the range of masses we have in our observed CMD is rather 
small and centered around 1 Mq . where the IMF is not especially sensitive to 
changes, we do not expect the effective SFH solution to significantly change 
if we assume another IMF. However, since the large number of lower mass 
stars are not well constrained by the data, different assumptions of the IMF 
will affect the normalization of the SFH, i.e. the total mass of the system. 
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Fig. 4. — (F435W - F555W , F555W) CMDs of field Fl (left-hand panel) and F2 (right-hand panel) obtained using deconvolved images. These contain 58143 
and 27963 stars respectively, and are calibrated onto the VEGAmag HST system. Note the difference between the CMDs in the region highlighted with dark 
blue boxes. The larger number of stars in Fl indicate the presence of a more significant intermediate-age population in this field compared to F2. The region in 
the box is not an actual "bundle" used in the derivation of the SFH but a similar one is used to obtain most of the information about the SFH of both fields: see 
Section|3]for more details. The light blue line indicates the 50% completeness level of our data in each field and the photometric errors from ASTs refer to a 
color ofTF435W-F555W) = 1. 



and metallicity bins in which the model CMD is to be divided. 
These simple populations represent the bins in which the SFH 
is to be determined. The boundaries of the bins that we used 
are [0, 0.5, 1, 2, 5, 14] x 10'' years in age and [0.02, 0.40, 0.80, 
1.00, 2.00, 4.00] X 10-2 in Z, corresponding to [M/H] w [-2.0, 
-0.67, -0.35, -0.25, 0.02, 0.32], assuming Zq = 0.019. These 
constitute 5 x 5 = 25 simple populations. The resolution in 
age and metallicity was selected, after several experiments on 
mock stellar populations, as the optimal choice for our data 
given the observational uncertainties'*. Note that the bin width 

Experiments on mock stellar populations were conducted as follows: we 
generated synthetic CMDs from arbitrary SFHs to which we apply the ob- 
servational effects of our observed CMDs obtained from the ASTs. We then 
solved for their SFHs using the lAC-POP/MinnlAC method as if they were 



in age increases significantly for older populations. This is 
due to the limits imposed by the crowding; we cannot extract 
more detailed information about the oldest stars. 

We then define five "bundles", macro-regions of the CMDs 
used for the fitting. We show the CMDs of Fl and F2 with 
the selected bundles superimposed in Figure l6] The bundles 
are subdivided into boxes, whose sizes vary from bundle to 
bundle. The bundles and boxes are equally sampled in the ob- 
served and model CMDs. Since the number of stars in each 
box is the information provided to the lAC-pop code, the dif- 

real data, adopting different age and metallicity resolutions as first reasonable 
guesses according to the limiting magnitude of our data. The SFH solutions 
were compared with the input SFHs, and when these agreed to within 1 — cr, 
we assumed the corresponding age and metallicity resolutions were optimal. 
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Fig. 5. — Left-hand panel: Hess representation with a logarithmic stretch of the synthetic CMD generated using lAC-STAR for a range of age between and 
14 Gyr and metallicities uniformly distributed at all ages between 0.0001 and 0.04. It contains 5 X 10* star's. Right-hand panel: Hess representation of the model 
CMD, i.e. the synthetic CMD after the observational effects have been simulated. It contains ~ 2 X 10* stars. The color bar indicates the number of stars per 
color-magnitude bin in logarithmic scale. This model CMD is the one to be compared with the observed CMD to derive the SFH. 
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Fig. 6. — Left-hand panel: CMD of field Fl in absolute magnitudes, assuming a distance = 24.53 and E(B-V) = 0.08, with the location of the bundles 
superimposed. Right-hand panel: As in the left-hand panel but for F2, assuming a distance fiQ = 24.45 and E(B-V) = 0.08. Each bundle in both CMDs is 
subdivided into boxes with sizes that vary from one bundle to another (see Tablel2|. This allows each CMD region used for the analysis to have different weights 
on the extracted SFH. Note that no bundles were added below the 50% completeness level of each CMD. Also note that most of the RGB and RC regions are 
likewise excluded of our SFH analysis: uncertainties in the physics governing these evolutionary phases are larger than those in the MS and SGB region. 
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TABLE 2 
CMD REGIONS USED FOR THE FITTING 



Bundle 


# of boxes 


Size of boxes (color, mag) 


CMD region sampled 


1 


500 


(0.01,0.20) 


lower MS 


2 


150 


(0.03, 0.30) 


upper MS 


3 


3 


(0.50, 0.40) 


SGB 


4 


7 


(0.50, 0.50) 


left of the MS 


5 


5 


(0.50, 0.90) 


Right of the RGB 



ferent bundle subdivisions provide the weights that a given 
CMD region has on the derived SFH. For instance, CMD re- 
gions well-populated and/or where the input physics is better 
understood (e.g., bundle 1) have smaller boxes than CMD re- 
gions where either the number of stars is smaller or the un- 
certainties in the input physics significantly impact stellar in- 
terior models (e.g., bundle 5). The properties of the boxes for 
each bundle are specified in Tablel2]and both the observed and 
model CMDs with one sample of the boxes superimposed are 
shown in Figure |7] Note that the boxes inside the bundles are 
shifted during the dithering process, as explained below, and 
only stars inside a bundle are considered in the analysis, no 
matter how big the box is. Also note that only stars brighter 
than the 50% completeness level were used to extract the SFH 
(cf. Fig.l6]l. Below this region, most of the information is lost 
and results obtained from lower-completeness regions are un- 
reliable (see also the right panel of Fig. [5]). Also, as mentioned 
above, we did not use most of the RGB and RC. Adding bun- 
dles in those regions significantly increased from ^ 2 to 
^ 5. Bundles 4 (bluer than the observed MS) and 5 (redder 
than the observed RGB) were adopted to mainly constrain the 
lowest and highest metallicity, respectively, in the observed 
CMD. There are nearly no observed stars in these regions (see 
Fig. |6]), whereas there are stars in the model CMD (see right 
panelof Fig.|5jl. The fact that stars of certain ages and metal- 
licities appear in the model CMD but not in the observed one, 
indicates that those simple populations are not present in M32. 

3) Solution. For a given parametrization, i.e., box sizes and 
simple population boundaries, MinnlAC counts the stars in 
each of the boxes for both the observed and model CMDs. 
The number of stars in each box is the input information to 
run lAC-pop. lAC-pop compares the observed and model 
star counts in e ach box using a modified merit-function 
(Mighell 1999 1, calculating which combination of simple 
populations best reproduces the observed CMD. A SFH so- 
lution is obtained as a linear combination of the simple pop- 
ulations. Thus, lAC-pop solves the SFH considering the age 
and metallicity as independent variables. 

4) Uncertainties and stability of the solution. To minimize 
biases in the solution due to the sampling of the CMD, Min- 
nlAC allows slight changes in the input parameters. The sim- 
ple populations (i.e. the age and metallicity bins) are shifted 
three times a 30% of their corresponding bin sizes for each 
of the following four different configurations: (1) shifting the 
age bin toward increasing age at fixed metallicity, (2) shift- 
ing the metallicity bin toward increasing metallicity at fixed 
age, (3) shifting both age and metallicity bins toward increas- 
ing values, and (4) shifting both age and metallicity bins to- 
ward decreasing age and increasing metallicity. Furthermore, 
for each of these 12 shifts, the boxes are shifted a fraction of 
their [color, magnitude] sizes three times: [80%, 20%], [20%, 
80%] and [20%, 0%], respectively. These 36 sets of param- 
eters are used to generate 36 individual solutions. The final 
SFH solution is the average of these. This "dithering" pro- 



cess significantly reduces fluct uations in the solution associ- 
ated with the sampling (Hidal go et al.|[20TT l. The standard 
deviation of the "dithers" provides a measurement of the un- 
certainties on the solution (see Aparicio & Hidalgo 2009 , for 
further discussion of uncertainties in the solution). 

To account for uncertainties in the distance modulus 
(±0.14: Paper I), reddening (±0.03: Burstein & Heiles 1982)), 
aperture corrections (Paper I), and other systematics possibly 
affecting the zero points of our photometry, we allow the ob- 
served CMD (not the model) to shift in both color and mag- 
nitude. The observed CMD is shifted four times in magni- 
tude and six times in color The bundles are correspond- 
ingly shifted. MinnlAC repeats the entire process of gen- 
erating the input information and averages the 36 individual 
solutions generated by lAC-pop, for each of the positions in 
a magnitude-color grid. The grid has 35 nodes, where the 
shifts in magnitude are (-0.14, -0.07, 0, 0.07, 0.14), and the 
shifts in color ai-e (-0.12, -0.09, -0.06, -0.03, 0, 0.03, 0.06). 
In total we generate 36 x 35 = 1260 individual solutions for 
each field (Fl and F2) and library (BaSTI and Padova/Girardi) 
combination. 

5) Final best solution. After the observed CMD-shifting 
and "dithering" process, we have 35 averaged solutions, one 
for each color-magnitude node. Among the 35 mean solu- 
tions, the one with the lowest is chosen to be the final 
solution that best reproduces our observed CMD. 

3.2. Best SFH solutions for Fl and F2 

As previously mentioned, for each shift in color and mag- 
nitude of the observed CMD, we average the 36 individual 
solutions as well as its corresponding xi- For a 35% binary 
fraction, the nodes at which the mean minimum xl,, i-e. xt min 
was reached were found at (6{F435W -F555W)[),6MfS55w) = 
(-0.09,0.07) for Fl with xl,,nin = 2-03, and at (6(F435W - 
F555W)o,6Mf555w) = (-0.06,0.0) for F2 with xl.min = 2-23. 
We consider the averaged solution corresponding to xl min 
the one that best reproduces our observations. Figure[8]shows 
how the mean xt varies as a function of the color and mag- 
nitude shifts applied to the observed CMD in Fl (left panel) 
and F2 (right panel), using the BaSTI library and 35% of bi- 
nary fraction. The crosses indicate the 35 nodes at which we 
calculated a mean x^, average of the 36 xi from the individ- 
ual solutions. The contours around the minimum xi (whose 
value is indicated in the figure at the position where it was 
found) show the 1-, 2-, 3- and 6-a confidence regions, with 
<T defined as the standard deviation of the mean ximin- We 
emphasize here that the shifts in the observed CMD at which 
we obtained the best solution do not necessarily represent cor- 
rections to the distance or reddening estimates, since photo- 
metric corrections and model systematics are also present. 

The xi min values suggest that the BaSTI library fits the data 
better than the Padova/Girardi isochrones for both fields (see 
Tablefflin the Appendix). Nevertheless, the solutions obtained 
with both libraries are very similar, with the Padova/Girardi 
isochrones generating a best-fit mean solution slightly more 
metal-rich than BaSTI. For simplicity, we consider the solu- 
tions obtained using the BaSTI library for most of the follow- 
ing analysis. 

4. RESULTS OF THE SFH ANALYSIS 
4.1. The SFH ofFl and F2 



The Star Formation History of M32 



9 



-2 



-1 




mm m 
Mwm liiii" 

i^iiii iiiiiiiiLi 

m in •«•' ••• 'WMi 

m i:: :: .. , i 



-0.5 



0.5 



1.5 



0.5 



1.5 



F435W-F555W 

Fig. 7. — CMD of field Fl (left-hand panel) and model CMD (right-hand panel) with the actual bundles and boxes used for the fitting superimposed. We show 
here one sample of boxes in which the stars are counted and compared between the observed and model CMDs. Note that the boxes inside the bundles are shifted 
during the dithering process, as explained in the text. Also note that only stars inside a bundle are considered in the analysis, no matter how big the box is. 

as in Fl. 

The most striking feature of Figure |9] is the significant star 
formation in Fl that occurred 2-5 Gyr ago. F2 is predomi- 
nantly old, with some contribution of young and intermediate- 
age stars from 0.5 to 5 Gyr ago, but its 2-5 Gyr old population 
is clearly not as prominent as that of Fl. We emphasize here 
that differences in the intermediate-age population between 
the fields were expected (see Paper 1 and Fig. |4]). However, 
the significant SFR in the 2-5 Gyr bin in Fl compared with 
F2 is rather surprising. As Fl has contributions from both 
M32 and M31 stars and F2 is expected to have a negligible 
contribution from M32, the derived SFHs suggest that the 2- 
5 Gyr old population in Fl is associated almost entirely with 
M32. We discuss this further in the next section. To obtain 
the mean age and metallicity of the system, we weight such 
quantities by the mass of each bin in age or metallicity, re- 
spectively. We call them hereafter mass-weighted mean age 
and mass-weighted mean metallicity. 

Figures 10 and 1 1 display the main results projected from 
the extractea SFHs of Fl and F2, respectively. We find that: 




-0.05 

6(F435W-F555W| 




-0.05 

6(F435W-F555WI 



Fig. 8. — Left-hand panel: Grid of color and magnitude shifts applied to 
the observed CMD of field Fl . The color-magnitude shift nodes are indicated 
with crosses. For each of these nodes, 36 individual solutions were obtained, 
and the average of its corresponding xi was calculated. We consider the so- 
lution at which the minimum mean xi, is reached as the best representation of 
our data. The xj, „j„ value obtained using BaSTl library for a binary fraction 
of 35% is indicated in the figure, at the position where it was found. Con- 
tours around this position show the 1-, 2-, 3- and 6-ct confidence regions, 
where a is defined as the standard deviations of the 36 xi- individual solu- 
tions. Right-hand panel: As in the left-hand panel but for F2. Note that the 
1 — (J confidence area is smaller in this case. 



In Figure|9]we show the best- fit mean SFH= 5'(f,Z) solution 
for Fl (top panel) and F2 (bottom panel) in a 3D-histogram 
representation, as well as the two projections ^'(f) (red line) 
and ^'(Z) (blue line), ^(f) is the SFR as a function of time or 
age distribution, i.e. 5'(f,Z) integrated over metallicity, and 
5'(Z) is the metallicity distribution function, i.e. 5'(f,Z) in- 
tegrated over time. Both distributions are normalized by the 
area in pc^. Recall that field F2 has ~ 1/2 the number of stars 



• Fl acquired 75% of its stellar mass between 5 and 14 
Gyr ago. Stars with ages of 2-5 Gyr contribute 23% 
of the mass in Fl. The remaining 2% of mass in Fl is 
constituted by stars younger than 2 Gyr. 

• Fl is metal -rich with an almost constant age- 
metallicity relation, to the limits of the age resolution 
of the CMD. 
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Fig. 9. — SFH(t,Z)= ^(t,Z) of Fl (top panel) and F2 (bottom panel) obtained using BaSTl models and assuming a 35% binary fraction. The blue and red lines 
are the two SFH projections: metallicity distribution ^(Z) and age distribution \I' (?), res pecti vely. Note that 'I'(Z) does not represent metallicity evolution, as it 
is integrated over age, and thus should not be compared with panel (d) of Figures [To] and pT] which show Z as a function of age. Each solution is calculated by 



aging the 36 solutions at the 



(5color grid (Sec. 



3.2 



- 2.03 for Fl and 2.23 for F2. Recall that the number of stars in F2 is ~ 1/2 of 
that in Fl . Note the prominent stellar population with ages 2-5 Gyr present in Fl but nearly absent in F2. Although differences were expected (note the different 
number of stars inside the blue box in Figure|4]and the results in Paper 1), the significant different SFRs in the 2-5 Gyr bin between the two fields is striking. 
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Fig. 10. — The SFH of Fl. (a) SFR as a function of time; (b) cumulative mass-weighted age distribution; (c) mass as a function of metallicity; (d) age- 
metallicity relation; and (e) comparison between the observed, calculated CMDs and a Hess representation of the residuals. The calculated CMD is obtained 
by randomly extracting stars from the model CMD in such a way that the final star distribution represents the calculated SFH. Both the observed and calculated 
CMDs were divided into the same 200 x 200 bins. The Hess diagram of the residuals in panel (e) shows the subtraction of the observed Hess diagram from the 
calculated one, in units of the Poisson uncertainties. The vertical solid line in panel (b) represents the mean age (~ 8 Gyr) of the system, and the dashed lines 
indicate the Icr deviation of that value. 
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TABLE 3 

Integrated quantities derived from the SFHs 



Field 



(Age) (Gyr) ([M/H]) (dex) int(SFH) (10* Mq) 



BaSTl library 

Fl 
F2 

M32(F1-F2) 
F2'' 



7.95 
9.12 
6.80 
9.15 



1.35 
0.80 
1.50 
1.27 



Padova/Girardi 
Fl 7.99 ± 1.33 

F2 9.03 ± 0.85 

M32(F1-F2) 7.03 ± 1.50 



-0.07 ± 0.10 
-0.19 ± 0.10 
-0.01 ± 0.08 
-0.10 ± 0.10 



0.01 ± 0.10 
-0.07 ± 0.10 
0.06 ±0.10 



5.17 ± 0.50 

2.59 ± 0.24 

2.60 ± 0.50 
2.50 ± 0.18 



5.88 ± 0.76 
2.81 ± 0.29 
3.07 ± 0.75 



^ SFH of F2 was derived u,sing BaSTI library with an extra age bin, from 5-8 Gyr. 

• Fl's mass-weighted mean age is 7.95 ± 1.35 Gyr and 
its mass-weighted mean metallicity is [M/H] = -0.07 ± 
0.10 dex (Table[3]). 

• F2 is predominantly old, with 95% of its mass already 
formed 5-14 Gyr ago. There is a small contribution of 
mass to the system after that, and it stopped forming 
stars ^ 0.5 Gyr ago. 

• F2 is also quite metal-rich, but is marginally more 
metal-poor than Fl, with a slight age-metallicity rela- 
tion showing a small increase in metallicity at younger 
ages. 

• F2's mass-weighted mean age is 9. 12 ±0.80 Gyr and 
its mass-weighted mean metallicity is [M/H] =-0.19± 
0.10 dex (Table[3]). 

The integrated quantities derived for the SFHs of Fl and F2 
using Padova/Girardi Library are also indicated in Table [3] 



Figures 10; and [TTt show comparisons between the ob- 
served (leftjand calculated CMD (middle) as well as the Hess 
diagram of the residuals in units of the Poisson uncertainties 
(right), for Fl and F2 respectively. The calculated CMDs have 
been obtained by randomly extracting stars from the synthetic 
CMDs in such a way that the resulting star distribution fol- 
lows the best calculated SFHs. For both Fl and F2, the model 
CMD shows reasonable agreement with the observed CMD 
throughout most evolutionary phases, which is also reflected 
in the residual Hess diagrams. The RC regions, however, 
show significant discrepancies. This is not surprising; due 
to uncertainties in, e.g., the mass loss during the RGB or the 
He content of the stars, that particul ar evolutionary stage is 
not well-modeled ( |Gallart et al.|2005"] l — ^but we have not used 
this region in deriving the solutions. There is also some dis- 
crepancies for magnitudes fainter than the 50% completeness 
level, but this region was also not used for the derivation of 
the SFHs. 

4.2. The SFH ofM32 as revealed by the lAC method 

To calculate the SFH of M32, we make use of the derived 
SFHs of Fl and F2\ Given flie fact that both SFHs have been 

^ We would ideally need a deep CMD composed solely of M32 stars to 
derive the SFH of M32, which we attempted to deiive in Paper I. Under the 
assumption that the M31 stellar populations in Fl and F2 are statistically the 
same, we subtracted the stars of the F2 CMD from the CMD of Fl taking into 
account the difference in crowding of the fields. This produced the deepest 
CMD of M32 yet obtained. However, the use of such CMD to extract the SFH 
of M32 would introduce uncertainties associated with the decontamination 
process. 



obtained using the same stellar population sampling, and as- 
suming that the SFH of M31 in Fl and in F2 is identical and 
that M32 is not present in F2, calculating the SFH of M32 is 
straightforward: we simply subtract the SFH of F2 from that 
ofFl. 



Figure 12 shows the inferred SFH of M32 for the first time 
calculatedTrom its resolved stellar population. We used the Fl 
and F2 SFHs shown in Figure |9] inferred using the BaSTl li- 
brary and a 35% binary fraction. We can see that a major burst 
of star formation occurred in M32 2-5 Gyr ago, responsible 
for ^ 40% of M32's current mass at Fl's location. This can be 
seen from the cumulative mass function, shown in panel (b) 
of Figure 13 Stars older than 5 Gyr contribute ~'55% of the 
total mass of M32 in this field. From this CMD-fitting anal- 
ysis, however, due to the limitations imposed by the crowd- 
ing of our fields, we cannot specify when the star formation 
started, whether it was constant over the 5-14 Gyr period, or 
if it peaked at some age. Integrated quantities derived from 
the calculated M32 SFH are indicated in Table [3] Note that 
the estimated mean age and metallicity of M32, ^ 6.8 Gyr 
and ^ -0.01 dex, respectively, are younger and more metal- 
rich that the mean age and metallicity of Fl, because M3rs 
mean age and metallicity in F2 is older and more metal-poor 
than M32 inFl. 

The age-metallicity relation for M32 is nearly constant 
(Fig. 131) although there is apparently a mild increase at ~ 
5 Gyrrollowed by a small decrease at ^ 2 Gyr. We note that 
an almost constant age-metallicity relation appears to suggest 
that M32 has not experienced any metal enrichment. How- 
ever, the lack of resolution in age means that we cannot ex- 
tract detailed information on stars older than 5 Gyr. Most 
of the chemical evolution of the system has likely occurred 
during that 5-14 Gyr period. M32's mass-weighted peak in 
metallicity is at [M/H] - 



:). 



'0.2dex(FigJ13; 
We show in the top panel of Figure 14j tne calculated CMD 



of M32, with its stars color-coded accormng to age. The CMD 
was obtained by randomly extracting stars from the model 
CMD, in such a manner that their star distribution follows the 
calculated SFH. This figure provides explicit information on 
the age interval that populates each region of the CMD as well 
as showing how the various ages combine. We see, for exam- 
ple, that stars of different ages contribute to the RC. Younger 
stars populate the brighter bluer portion of the RC while older 
stars populate the fainter, redder portion of the RC. The BP is 
only populated by stars younger than 2 Gyr. The bottom panel 
shows the CMDs produced by each age interval considered in 
the extraction of the SFH. We can appreciate in detail the dif- 
ferences between each CMD as the ages vary, from only an 
extended main sequence (bottom left panel, ages ^ 0.5 Gyr) 
to a CMD with well-populated RGB, RC and AGB evolution- 
ary phases (bottom right panel, ages of 5-14 Gyr). Note the 
presence of only few BHB stars in the bottom right panel, as 
expected for systems as metal-rich as M32; in the composite 
CMD (top panel), these few BHB stars are mixed with young, 
blue stars in the extended MS. 

Finally, we can qualitatively compare these results with the 
SFHs derived for some of the Local Group dwarf satellites, 
which have been analyzed by the same or very similar meth- 
ods'". It is interesting to note that the mean age derived for 
M32 is comparable to that of dwarf irregular galaxies, s uch as 
the Small and Large Magellanic Clouds and Pegasus (Noel 



* A more extensive and quantitative comparison of our results with other 
the SFHs of Local Group satellites is beyond the scope of this paper. 
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Fig. 12. — SFH of M32 obtained after subtracting the calculated SFH of F2 from that of Fl. We find two dominant populations contributing to the SFH of 
n52. One IS 2-5 Gyr old and contributes ~ 40% of the total mass of 1V132 at Fl. The population older than 5 Gyr contributes ~ 55% of the total M32's mass at 
V\ . Note that some of the stars y ounger than 2 Gyr are quite metal-poor compared to the nearly solar mean metallicity of M32. This suggests that these are BS5 
and may be the first direct evidence of such a population m M32. 



[eTaL] [20091 |Hams & Zan ts^lOOOl |Dolphin et aL][2005l l, 
whereas it is quite young compared with the mean ages of 
dwarf spheroidal galaxies, such as Cetus, Tucana, and Ursa 



Minor ( [Monelli etal 1201 Oa|bl [Dolphin et al.|2005| l. We also 
note that a synthetic CMP analysis was performed on archival 



WFPC2 data of M32 by 'Dolp hin eFaLl ( |200? ). Their results 
suggest an older mean mass-weighted age for M32, ~' 8.5 
Gyr However, the WFPC2 data not only are shallower than 
our data but also contain significant contamination by M31 
stars, which were not taken into account in their SFH analy- 
sis. 

4.2.1. Exploring the SFH solution and its robustness 

We now address the robustness and uniqueness of the SFHs 
derived here. First, due to the complex parameter space and 
thus multiple local minima involved in the process to find a 
SFH, an algorithm that guarantees finding the global mini- 
mum is strongly desirable. As mentioned above, we have used 
a genetic algorithm to find the minimum (see Charbonneau 



[1995| l. This type of algorithm, unlike the so-called 'down- 
hill' algorithms, is designed such that the solution found is in- 
finitesimally close to or at the global minimum independently 
of the initial seed that started the process, provided that a suffi- 
ciently large number of generation s (i.e. variations of individ- 
uals and mutations) is performed. jAparicio & Hidalgo[ ( |2009| l 
have tested this and found that w 10^ generations are enough 
to guarantee that the solution will be at or as close as possible 
to the global minimum. We have performed 2x10^ gener- 
ations per solution, which assures us that we have reached 
convergence. 

In addition, the lAC method does not introduce any sys- 
tematic error to the SFH solution, provided that the age and 
metallicity bins used to extract the SFH are appropriate to 
the observed CMD. This has been verified by several tests on 
mock stellar populations performed at Instituto de Astrofisica 
de Canarias in which the SFH of mock galaxie s of known 
SFH have been recovered rather accurately (see [Aparicio & 
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Fig. 13 . — The SFH of M32. (a) SFR as a function of time, clearly indicating the two dominant populations: at ~ 8 Gyr and ~ 4 Gyr; (b) cumulative mass- 
weighted age distribution which shows how much each population contributes to the total mass of M32 at F2; (c) mass as a function of metallicity. indicates tEe 
mean metallicity of the system, roughly solar; and (d) age-metallicity relabor^ieariyjxmslary^^l^ 



oF^^TTn Fl. I'he dashed lines indicate the Icr deviation of tins value. 



|Hidalgo|[2009| [Hidalgo et d:][2009l ^Un) . In this work, we 
have also performed several experiments on mock data to find 
the appropriate age and metallicity resolution at which, ac- 
cording to the observational effects of our observed CMDs, 
the SFH solutions of mock galaxies are within l-a to the 
input ones (see Sec. 3.1 above). Thus, each SFH solution 
obtained is "unique," by which we mean that combinations 
of simple populations within the error bars of the SFH will 
produce CMDs indistinguishable from the best fit CMD. Any 
other SFH which is combination of simple populations signif- 
icantly different that those of the final SFH (i.e., not possible 
within the error bars of our solution) will produce a CMD sig- 
nificantly different than the best-fit CMD and, therefore, than 
the observed one. 

Finally, even though we only show our best solution, we 
have explored other solutions that give similar good fits. We 
found that the SFH of the solutions at l-cr confidence area 
(see Fig.|8]l are not significantly different than the best one, for 
both Fl and F2. The mass percentages per age range slightly 
vary from one solution to the other, but the overall SFH re- 
mains the same. Taking these nearby solutions into account, 
we find that M32 at Fl has 40% ± 17% of its mass in a 2-5 
Gyr old, metal -rich population and ^ 55% ± 21% of its mass 
in stars older than 5 Gyr, with slightly subsolar metallicities. 
The uncertainties represent the l-cr error of our best solution 



and the variations of these percentages when considering so- 
lutions of similar good quality fit. 

4.2.2. Young population (Ages < 2 Gyr) vs. Blue Stragglers 

In Paper 1, we discussed the possibility that the fainter stars 
in the BP of M32 could be old BSS rather than a young stellar 
population with ages < 2 Gyr. However, the analysis pre- 
sented in Paper 1 did not allow us to confirm or rule out either 
case. BSS are stars hotter, bluer and brighter than the MSTOs 
in a CMD, thus generating a blue plume. Given their loca- 
tions on the CMD, they are burning hydrogen in their cores 
with masses larger than the turn-off mass, which indicates 
that some sort of mechanism rejuvenated their inner layers. 
Although such a mechanism is still a matter of debate, there 
are currently two theoretical possible scenarios to explain the 
BSS origin: they are the r esult o f either a collision between 
stars (e.g., Sigurdsson et a l.|1994 1 or mass-trans fer in a binary 
system (e.g., ,McCrea|1964irarney et al.|2001[ | 



We investigate the nature of these stars from the SFH pre- 
sented here. Stars younger than 2 Gyr constitute ^ 4% of 
the total mass of M32 at Fl. Figure [12] shows that some of 
the young stars, produced by a very low SFR event at look- 
back times < 2 Gyr, are rather metal poor ([M/H] ~ -0.7) in 
comparison with the mean metallicity of M32 ([M/H] ~ 0.0). 
Given the almost constant age-metallicity relation for M32 
and the presence of intermediate-age stars (2-5 Gyr old) of 
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Fig. 14. — Top panel: Calculated CMD of M32, obtained by randomly extracting stars from the model CMD in such a way that they follow the derived SFH of 
M32. The stars are color-coded according to age, as indicated in the bottom panel, except for 2-5 Gyr old and stars older than 5 Gyr, shown as yellow and black 
dots, respectively, in the top panel whereas gray-scale Hess representations of their CMDs are shown in the bottom panel. Note how the various ages fit together 
and the age interval that populates each region of the CMD. Bottom panel: Each CMD is composed by stars of a different age interval, from only an extended 
main sequence (left panel, ages ~ 0.5 Gyr) to a CMD with well populated RGB, RC and AGB evolutionary phases (right panel, ages of 5-14 Gyr). Note the 
differences in the MSTO region and fainter MS in the last two CMDs. The MSTOs for the younger population (2-5 Gyr, yellow dots in the top panel and a Hess 
representation of their CMD in the bottom panel) are brighter and bluer than the ones for the 5—14 Gyr old population (black dots in the top panel and a Hess 
representation of their CMD in the bottom panel). 
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solar or even higher metallicity, it is unlikely that M32 con- 
tains at the same time younger stars with significantly sub- 
solar metallicities. The most plausible explanation is that 
these stars are BSS belonging to an old metal-poor po pula- 
tion. BSS are found in o pen and globular clusters ( Ferraro 
eraDl2004; Mapelli et al. 2004, 2006, Piott o et al. 2004; de 
M'archi et al. 2006), dwarf spheroid aTgalaxies (Hurley-Keller 
et al. . l 999; Carrera et al. 2002; M omany et al.|2007^ . Mapelli| 
et al.|2009; MonelH et al. 2010a), and even in t he Milky Way 



halo field population ( |Preston & Sneden|2000 V Therefore, it 
seems natural to consider that they can also be found in an 
elliptical galaxy. These stars represent ^ 2% of the mass of 
M32 in Fl and might be the first direct evidence of BSS in this 
galaxy. An alternative explanation could be that these young 
and metal-poor stars were generated by an episode of late in- 
fall of metal -poor gas. However, if we assume that M32 is 
interacting with M31, we would not expect M32 to accrete 
gas, but instead to lose gas to M3 1 through stripping. 

The other ^ 2% of stars with ages < 2 Gyr that we find 
in the SFH inferred for M32 may indeed represent a young 
metal-rich population in M32 at Fl. 

5. THE STAR FORMATION HISTORY OF M32 

By combining the results in the present work with the anal- 
ysis in Paper 1, we can finally provide a detailed and complete 
SFH of M32. We conclude that M32 has had an extended 
SFH and is composed of two main dominant populations at 
Fl: 40% ± 17% of the mass in a 2-5 Gyr old, metal-rich 
population and ^ 55% ± 21% of the mass in stars older than 
5 Gyr, with slightly subsolar metallicities. We confirm the ex- 
istence of the younger (< 5 Gyr) stars through the presence 
of bright AGB stars observed in Paper 1, with the appropriate 
ages. From the RC, RGB bump and AGB bump analyzed in 
Paper I, the bulk of the old population is 8-10 Gyr old. We 
therefore do not expect a significant contribution from stars 
older than 10 Gyr in M32 at Fl. Nevertheless, there is a hint 
of the presence of a few ancient metal-poor stars present in 
M32, as revealed by the 2-a detection of RR Lyrae belonging 
to M32 at Fl. The remaining ^ 4% of the mass is roughly 
equally divided between a young metal-rich population and 
a young metal-poor population. We associate the latter with 
blue straggler stars belonging to an old (likely metal-poor) 
population. 

The age-metallicity relation for M32 is nearly constant 
within our age resolution, although there is a small increase 
in metallicity at at ^ 5 Gyr followed by a small decrease 
at ages younger than ^ 2 Gyr The mass-weighted mean 
metallicity of M32 is [M/H] -0.01 dex with a peak at 
[M/H] ^ 0.02dex. We emphasize here again that an almost 
constant age-metallicity relation appears to suggest that M32 
has not experienced metal enrichment; but as in Fl, this is due 
to the poor age resolution and does not imply the lack of an 
age-metallicity relation. Stars with metallicities lower than 
[M/H] < -Idex only contribute ^ 5% of the total mass of 
M32 at ^ 2' from its center. This is consistent with the pho- 
tometric metallicity function (MDF) of M32 derived in Paper 
I, which shows that the majority of the stars has a slightly 
sub-solar metallicity at [M/H] ^ -0.2 dex. The MDF also in- 
dicated that metal-poor stars with [M/H] < -1.2 contribute 
very little, at most 6% of the total V-Ught to M32 or 4.5% 
of the total mass in Fl, implying that the enrichment process 
largely avoided the metal poor stage. 

5.1. On the integrated light ofM32 



The results obtained in this work are a fundamental step 
for understanding the formation and evolution of other low- 
luminosity spheroidal systems (elliptical galaxies or bulges). 
Since integrated light spectra are, in general, the only means 
available to study the stellar populations of these galaxies, we 
strongly rely on unresolved stellar population models to learn 
about their SFHs. These models, which have become very so- 
phisticated in disentangling the non-trivial age and meta llicity 
degeneracy (Worthey 1994; Bruzual & Chariot 2003; Thomas 
|et al.|200^ [Vazdekis et al.|2010 ), still suffer from several un- 
certainties: e.g., it is difficult to distinguish between a young 
or hot old population since the latter is not necessarily ac - 
counted for in the models (e.g. Maraston & Thomas 2000). 
Calibration of these models, which requires observations of 
individual stars in elliptical galaxies, is a key ingredient that 
needs to be further developed. As briefly mentioned in Sec- 
tion 1, M32 is located at a distance such that both integrated 
spectroscopy and photometry of its individual stars can be 
studied in great detail. A comparison of the stellar parame- 
ters obtained using resolved stars and integrated luminosity is 
fundamental to provide a calibration to the unresolved stellar 
models with an actual elliptical. 

Extensive spectroscopic studies of M32 have been per- 
formed, mostly in the central reg i ons and out to ^ 1 (e.g., 
[19801 
1 studi 



Gonzalez 



T993j |Worthey|[2004l [RoseeFa: 



O'Connell 

2005p . All studies agree that tlie central stellar population 



has an SSP-equivalent age of 2.5-5 Gyr and roughly solar 
metallicity, with an age gradient that increases the age at 1 r^. 
by ^ 3 Gyr and a mild negative metallicity gradient. Vari- 
ous integrated-light studies have suggested that M32 under- 
went a period of significant star formation in the recent past. 



i.e., about 5- 8 Gyr ago, (e.g., |0'Connell| 1980[ |Pickles | 1985 
Bica et aL] l990) based on the presence of enhanced nj3 ab- 



sorption m the integrated spectrum of M32, a signature of 
an intermediate-age popula tion (e.g. ,|Rose|1994||Trager et al.| 
2000, Worthey 2004; Schi avon et at 2004; Rose et al. 2005[ 
,Coelho et al..200 9). To date, only Coelho et al. (2009) have 
attempted to probe the unresolved stellar populations as far 
from the center of M32 as the ACS/HRC field presented in 
this paper lies, using longslit observations with GMOS on 
Gemini. They propose that an ancient and intermediate-age 
population are both present in M32 and that the contribution 
from the intermediate-age population is larger at the nuclear 
region. They claim that a young population is present at all 
radii, and they further suggest that there is a strong compo- 
nent of either very young (< 0.3 Gyr) and/or very old (> 10 
Gyr), metal-poor stars even in their outermost field. 
We can use the inferred SFH of M32 to compute line 



strength indices using the models of Bruzual & Chariot ( 2003 
hereafter BC03) and to calculate single stellar population 
(SSP)-equivalent parameters that can then be compared with 
the values obtained from the integrated light of this galaxy. 

Using the BC03 models, we obtain a B-band luminosity- 
weighted mean age and metallicity of 4.9 Gyr and [M/H] = 
-0.1 2 dex, respectively, for M32 at Fl from its resolved SFH. 



Coelho et al.|((2009i) find an average luminosity-weighted age 
7± L5 GT^ 



of 5.7 ±1.5 Gyr using BC03, which agrees with our result 
within the uncertainties, but their inferred mean metallicity is 
much lower, [M/H] =-0.6 ±0.1 (see their Table 3). More- 
over, as mentioned above, they suggest that there is a strong 
component of either very young (< 0.3 Gyr) or very old (> 10 
Gyr), metal-poor stars in their field at a radius similar to our 
Fl field, which is inconsistent with our data. 
We have also calculated SSP-equivalent parameters that 
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an SSP-equivalent age of ~ 8.4 Gyr and a SSP-equivalent 
metallicity of [M/H] - -0.13 dex. Figure p^lshows the SSP- 
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Fig. 15. — SSP-equivalent age (top panel) and metallicity (bottom panel) 
values infeiTed from line-strength fits given by Worthey (20041 as a function 
of radius in M32. The lines are linear fits to the data. We can see the steep 
positive age gradient from M32's center to 1 ~ 40". An extrapolation of 
these fits to log(110") ~ 2.04 gives the SSP-equivalent age and inetallicity 
values at Fl : 8.9 Gyr and [M/H] = -0.23 ± 0.03 dex, respectively. The SSP- 
equivalent parameters calculated from the infeiTed SFH of M32 (blue dots) 
are in stark contrast with predictions of spectral studies. 

can be compared with the values obtained from integrated 
spectra of this galaxy. Using the BC03 models, the SSP- 
equivalent values of M32 obtained from its inferred SFH at Fl 
are 2.9 ± 0.2 Gyr and [M/H] = 0.02 ± 0.01 dex, respectively. 
Given the radial li ne-strength gradients present in M32 (e.g., 
[Rose et aL]|2005| l, we cannot directly compare these SSP- 
equivalent values with t hose obtained from central integrated 
spectra of M32 (by, e.g., Trager et al.|2000[ ). We therefore do 
the f ollowing. We use the values of the line strength indices 
from|Worthey|([2004 ) and compute the SSP-equivalent param- 
eters from polynomial fits to the absorption-line strengths as 
a function of radius (his Table 1) using the modified BC03 
models described in Trager et al. (2008 ). We then fit straight- 
line gradients to the SSP-equivalent parameters as a function 
of radius and extrapolate these fits to 110", Fl 's position. The 
SSP-equivalent age and metallicity of M32 at Fl from this 
extrapolation are 8.9 ±0.5 Gyr and [M/H] = -0.23 ±0.03 
dex, respectively. We note that Wort hey ' s values for Mgb ar e 
low compared with 'Gonzalez ( 1993| l and |Trager et al. ( 1998| l, 
and thus the SSP-equivalent age we have obtained may be 
slightly overestimated whereas the SSP-equivalent metallicity 
may be underestimated. Taken this into account, we obtained 



equivalent parameters fro m [Worthey ( |200f| 7The linear fits to 
the log(age/Gyr) and [M /H] as a function of log(r/") are also 
shown. The blue stars indicate the SSP-parameters obtained 
from the infeiTed SFH of M32 from BC03 models. 

The predicted age and metallicity at Fl from the extrap- 
olation of the absorption-line gradients are much older and 
more metal-poor, respectively, than those obtained from the 
inferred SFH of M32. This suggests that either the extrapola- 
tion of line strength indices or the stellar population models, 
or both, may be in error, but we are currently unable to discern 
which. Color profiles in many colors of M32 are rather flat 
( Peletier|1993 1. Since M32 does not contain dust, integrated 
colours can be good population indicators and the fact that 
there are no g radients in colors ag ree with the results from the 
inferred SFH. Davidge & Jensen (2007 ) have also challenged 
the radial gradients in mean stellar parameters obtained from 
spectral studies. They find no evidence for a radial age gra- 
dient in M32, based on the properties of observed brig htest 
AGB stars, in contrast to the results by Worthey ( |2004[ ) and 
|Rose et ari ('2005 ), who found (as described above) a signif- 
icant radial gradient in the mean luminosity-weighted age of 
M32. 

To farther investigate this apparent contradiction, in tegral- 
field spectroscopic observations with VIRUS-P (Hill et al. 
2008|l have been taken at Fl and F2 and will be analyzed in a 



future paper. This will provide spectra of the integrated stellar 
light of M32 for the fundamental calibration for the study of 
stellar populations. 

5.2. On the formation ofM32 

Certainly, the most striking result of this work is the sub- 
stantial contribution of 2-5 Gyr old metal-rich stars to the to- 
tal mass of M32 at Fl. How has an elliptical galaxy like M32 
formed such a young population of stars? What is the ori- 
gin of this population? In this section we attempt to address 
these questions and discuss, in particular, the most popular 

proposed formation scenarios for M32. 

A form ation scenario for M32 has been proposed by KorJ 
mendy et al. (2009 , hereafter K09), in which M32 is a normal7 
low-luminosity elliptical galaxy. K09 find that both central 
and global parameter correlations from recent accurate pho- 
tometry of galaxies in the Virgo cluster place M32 as a nor- 
mal, lo w-luminosity elliptical galaxy in all regards. K09 fit a 
Sersic|(1968 ) profile to the SB of M32 with n = 2.8, in agi-ee- 



ment with Sersic indices of other low-luminosity ellipticals 
studied by K09. They interpret the light at the center of M32 
that was not fit by their Ser sic profile as a signature of forma- 
tion in dissipative mergers ( Mihos & Hernquist|1994| l. Extra 
central light is a general feature of coreless galaxies and is 
observed in all the other low-luminosity ellipticals of K09's 
sample. Within this scenario, the metal-rich 2-5 Gyr old stars 
contributing to w 40% ± 17% of M32's mass at Fl could be 
the result of such a dissipative merger event. Thus, the pro- 
genitors of M32 should have been very gas-rich spiral galax- 
ies, like M33 for example. However, such progenitors should 
have stellar masses of the order of Iff" Mq, whereas M33's 
stellar mass is 3 x \(f Mq. There are in fact no gas-rich 
spiral galaxies near M31 of the appropriate stellar mass. 

An altern ative scenario for th e formation of M32 has been 
proposed by Bekki et al. pOOlj hereafter BOl), who assumed 
that M32 is the result of a low-luminosity spiral galaxy, whose 
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bulge, unlike most of its outer disk, survived its dynamical 
interactions with M31^. In their N-body/smoothed particle 
hydrodynamics simulations, BOl considered a gas-rich low- 
mass disk galaxy with a bulge orbiting a massive disk galaxy 
like M31. About 0.75 Gyr after the interactions have started, 
the outer stellar disk (from 2 kpc to 5 kpc) of the spiral galaxy 
is stripped away and only keeps « 40% of its initial mass in 
stars initially located in the central regions, i.e., within 2 kpc 
of the center New star formation is triggered by the interac- 
tion of the gas-rich spiral with M31 but the outer new stel- 
lar component is also tidally stripped away, and consequently 
only the central starburst component survives. On the other 
hand, the bulge is only weakly affected by tidal interactions 
with M31 due to its compactness, and only w 19% of its mass 
is lost. At the end of their simulations, there is a fractional 
disk, bulge, and new stars mass ratio of w 49%, ~ 42%, and 
w 0.9%, respectively, within 2 Kpc of the remnant compact 

galaxy. Our field Fl is located at 110", i.e. - 0.5 kpc from sentation of the new SFH solution for F2 



and B06 derived SFHs at each field in great detail. Differ- 
ences between these SFHs were mainly found in the age and 
metallicity distributions of stars older than 5 Gyr Within this 
age range (5-14 Gyr) we do not have the resolution required 
to inspect different bursts of star formation in F2 in detail, 
given the SFH extracted in Section!?] We can, nevertheless 
study the SFH of F2 in more detail tnan what is presented in 
Section |4] As we show in Figure ffl the CMD of F2 is 0.5 
mag deeper than the one of Fl, wmch allows us to obtain in- 
formation of fainter, i.e. older, MSTOs at F2'^. We therefore 
extracted again the SFH of F2 following the steps indicated in 
Sectionl3] but with an extra bin in the age, from 5 to 8 Gyr, for 
the simple populations considered. The boundaries of the bins 
in age are in this case [0, 0.5, 1, 2, 5, 8, 14] x 10^ years. The in- 
ferred best-fit mean SFH of F2 with this new resolution in age 
was found at {S{F435W -F555W)o,SMfS55w) = (-0.03,0.00) 



with ximin - 2.12. Figure 16 shows a 3D-histogram repre 



the galactic center at M32's distance and, assuming either an 
insi de-out or outside-in formatio n scenario for the disk (see 
e.g., Sommer-Larsen et al.|2003[ and discussion in Section[6]l 
and considering that we are looking at a ^ 0.5 Rd, where 
= 0.9 kpc is the scale length radius of BOl's disk, the disk 
stars that we should be observing there would have mean ages 
of 8-12 Gyr Assuming that the bulge is predominantly old, 
this scenario is difficult to reconcile with our results, given the 
substantial 2-5 Gyr old intermediate-age population detected 
in this work. However, there has been some indications of 
small bulges which could have extended SFHs, similar to that 
of M32, with 10-30% of their total mass at look-back times 



between 0.5 and 5 Gyr (Thomas & Davies 2008). 

Therefore it is unclear from our SFH results what the pre- 
ferred model for the origins of M32 is. While specific origin 
models differ in detail, the general outlines overlap enough to 
make choosing a specific model difficult with the age reso- 
lution of our current SFH. More observations of M32-analog 
systems and simulations of spheroidal systems with similar 
SFHs to that we have presented are needed to shed light on 
M32's origins. Furthermore, finding evidence of a stellar 
stream in the halo of M31 with the ages and metallicities ob- 
tained for M32, which should be left if a major stripping of 
M32 by M3 1 has occurred, would help to constrain the mod- 
els and assess the validity of the proposed scenarios. 

6. THE DISK AND SPHEROID POPULATION OF M3 1 IN F2 

In Paper I, we compared our findings in F2, in particular its 
metallicity distribution function (MDF), with several previ- 
ous works on the disk and bulge of M31 (e.g. [W illiams 2002; 
[Worthey et ai:p005|[Olsen et al.|2006HBrown et al..200 6 ). We 
found a reasonably good agreernent with most studies. In this 
section we discuss our new, quantitative results on the stellar 
populations at F2 and their implications for the formation of 
the M3rs disk. M31 seems to have formed most of its stars 
between 5-14 Gyr ago at F2. As mentioned above, we cannot 
precisely indicate when the star formation started in either Fl 
nor F 2 but we can see that M31 is older than M32 in Fl. 

Bro wn et aL] ( |2006| hereafter B06) analyzed three CMDs 
of different regions of M3 1 : the spheroid, stream and outer 
disk. These CMDs reached well below the oldest MSTOs, 

' The idea that M32, as well as other small high-surface brightness galax- 
ies, is a tidally truncated galax y has been discussed several decades before 
BOl models. In, for example, Faber]{T973|, the original truncated galaxy was 
a more massive elliptical galaxy, from which only the tightly bound core of 
the original elliptical remains after a strong tidal interaction. 



We can now dis- 
tinguish two main populations that contribute substantially to 
our background field F2, instead of only one old population: 
ss 30% ±7.5% of the total mass in F2 is composed of a 5-8 
Gyr old metal-rich population and w 65% ± 9% of the mass 
is composed of a 8-14 Gyr old, metal-poor population. An 
age-metallicity relation shows a slightly steeper slope from 
an old metal-poorer population to younger metal-richer ones 
than before, as shown in Figure [TT] We are still not able to an- 
swer when the star formation started in F2. Nevertheless, our 
results for the mean age and metallicity for F2, 9.12 ± 1.21 
Gyr and -0. 10 ± 0. 10 dex respectively, are in good agreement 
with B06 results for their outer disk field, 8.5 Gyr and -0.4 
dex, respectively''. In addition, young stars, with ages be- 
tween 0.3 and 1 Gyr, that populate the BP in the CMD of 
F2 do not contribute significantly to the total mass, which is 
also in agreement with B06's results. Interestingly, kinematic 
data in our field imply that both the disk and spheroid of M3 1 
contribute to the populations in F2 (K. Howley, 2010, priv. 
commm). This was also the case for the outer disk field of 
B06. B06, however, attempted to disentangle both popula- 
tions assuming that their spheroid field was representative of 
the spheroid population present in their outer disk field. By 
subtracting the spheroid population, they obtained a younger 
mean age for the outer disk of M31 — but still older than 5 
Gyr 

Given the resolution allowed by the depth of our data, the 
inner disk and spheroid populations of M31 (at 5 kpc from 
its center) seem to be indistinguishable from the outer disk 
and spheroid ones (at 25 kpc from M31 galactic center, B06). 
Even though we are unable to subtract the spheroid popula- 
tion that contributes to our field F2, most likely the mean age 
of M3rs disk at F2 is younger than 8.72 Gyr and older than 5 
Gyr, given the negligible contribution of stars younger than 5 
Gyr This result supports the inside-o ut disk formation model s 
by e.g., Abadi et al. (2003a b); So mmer-Larsen et al.| ( |2003[ ). 
Abadi et al. find a mean age of 8-10 at 2 kpc, which radially 
decreases to 6-8 at 20 kpc. Sommer-Larsen et al. simulated 
two spiral galaxies, with two different scenarios of disk for- 
mation: inside-out and outside-in. Our expected mean age 

The previous selection of age and metallicity bins to derive the SFHs 
was strictly based on the resolution imposed by the CMD of Fl. In order to 
subtract the SFH of F2 from that of Fl, we required the simple populations 
considered be exactly the same. 

' The cited values correspond to the results obtained by B06 when a 40% 
binary fraction was assumed. 
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Fig. 16. — A more detailed SFH of F2 in a 3d-histogram representation than that shown in Figurel9] This SFH of F2 was constructed this time with an extra bin 
in age covering 5-8 Gyr We now find two dominant populations of M31 at F2: An old more metaPpoor population, older than 8 Gyr, and an intermediate-age 
more metal-rich population, 5-8 Gyr old. Stars younger than 5 Gyr old only contribute ~ 5% of the mass of M31 at F2. 



for the disk of M31 at F2 agree with both scenarios within 
their unc ertainties, assuming a stellar disk scale length of ~ 5 



kpc (e.g., |Walterbos & Kennicutt|1988||Worthey et al.|2005b . 
They find that, at 1 disk scale length, the mean ages of both 
simulated disks are ^6-8 Gyr. However, the significant frac- 
tion of stars younger than 5 Gyr predicted by their outside-in 
model at F2 is not supported by our data. Thus, we favor 
their insi de-out model. Furthermore, th e inside-out formation 
model of Sommer-Larsen et al.| ( |2003| l predicts that the disk 
has almost no age gradient which, although surprising, is also 
in agreement with the comparison of our and B06 results at 
different disk locations. They explain that this prediction is 
a consequence of the non linear dependence of the SFR on 
the cold gas density, which makes the SFR rather low in the 
outer disk at late times, thus the average outer disk stellar age 
is quite high. An alternative scenario for the absence of an 
age gradient, found when comparing our results with those 
of B06, is the radial migration of stars seen in r ecent simula- 
tions of isolated disk form ation and evolution ( |Roskar et al.| 
|2008[ [Minchev et al.|2011 ). In these simulations, inside-out 
disk growth yields a negative age gradient within the break ra- 
dius (2-3 disk scale length), after which there is a positive age 
gradient due to the secular redistribution of stars, given their 
interactions with transient spiral density waves. Of course. 



we should keep in mind that what we presented here are the 
results of one field in the inner regions of M3 1 and we need 
more observations and statistics to either confirm or rule out 
what we suggest. The multi-cycle Panchromatic Hubble An- 
dromeda Treasury (PHAT) project, which will cover 1/3 of 
M31 with HST WFC3 and ACS observations, will resolve the 
SFH of the disk of M3 1 : our observations and analysis merely 
hint at what PHAT is likely to find. 

7. SUMMARY AND CONCLUSIONS 

We used deep HST ACS/HRC observations to derive the 
SFH of M32 for the first time from a detailed modeling of its 
CMD. The two fields observed, one closer to M32 (Fl) and a 
background M3 1 field (F2), were introduced and used in Pa- 
per I to construct deep CMDs of Fl and F2, and the deepest 
optical CMD of M32 yet obtained. The lAC-POP/MinnlAC 
method was used here to compare the distribution of stars in 
the observed CMDs of Fl and F2 with that of a model CMD. 
We obtained the SFH of M32 by linearly subtracting the SFHs 
of F2 from that of Fl . The use of different stellar evolutionary 
libraries (BaSTI and Padova/Girardi) and assumptions of bi- 
nary fractions (0, 0.35, 0.7, and 1) did not significantly modify 
the solutions obtained, indicating that our results are robust. 

Combining our present results with those of Paper I, we 
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Fig. 17. — The SFH of F2 infen'ed using an extra bin, i.e. more resolution, in age. (a) SFR as a function of time; (b) cumulative mass-weighted age distribution; 
(c) mass as a function of metallicity; (d) age-metallicity relation. The vertical solid and dashed lines in panel (b) represent the mean age (~ 9.15 Gyr) and Icr 
deviation of that value, respectively. 



provide an unprecedented census of the stellar content of 
M32. The derivation of the SFH presented in this paper is 
independent of the analysis performed in Paper 1. In spite of 
using the same data, the CMD regions that we have probed 
in this work are largely different from those used in Paper 1. 
Our analysis of these regions have allowed us to obtain de- 
tailed information about the young and intermediate-age pop- 
ulations of M32, whereas only the broadest sketch of these 
populations were possible in Paper 1. Conversely, detailed in- 
formation about the older populations cannot be obtained with 
our current approach, and therefore we rely on the qualitative 
results of Paper 1 for those populations. 

The main finding of this work is that M32 is composed 
of two main dominant populations at Fl: ^ 40% ± 17% 
of the mass in a 2-5 Gyr old, metal-rich population and 
~ 55% ±21% of the total mass in stars older than 5 Gyr, 
with slightly subsolar metallicities. Its mass-weighted mean 
age and metallicity ai-e (Age) = 6.8 ± 1.5 Gyr and ([M/H]) = 
-0. 1 ± 0.08 dex, respectively. Even though we are unable to 
specify when the star formation started in M32 at Fl, we make 
use of the analysis of Paper 1 to constrain the older population. 
We know from the RC, RGb and AGB bumps that the bulk of 
the old population is 8-10 Gyr old. Thus, we do not expect a 
significant contribution from stars older than 10 Gyr in M32. 
There has been, however, a marginal detection of RR Lyrae 
belonging to M32 at Fl, which reveal the presence of a few 



ancient metal-poor stars in M32 (|Fiore ntino et al.|2010 1. The 
remaining ^ 4% of the mass is distributed in geriuine young 
metal-rich stars 2%) and young metal-poor stars (~ 2%) 
which we associate with blue straggler stars belonging to an 
old metal-poor population. In addition, we used the inferred 
SFH of M32 to calculate the SSP-equivalent age and metal- 
licity parameters from unresolved stellar population models, 
which ai-e 2.89 ± 0.15 Gyr and [M/H] = 0.02 ± 0.01 dex, re- 
spectively. These values, however, contradict spectroscopic 
studies, which show a steep age gradient from M32's center 
to 1 r^. 

Based on our present results, it is not currently possible to 
choose a preferred model for M32's origins between two pop- 
ular ones: a true low-luminosity elliptical or a former spiral 
galaxy whose bulge survived its dynamical interaction with 
M31. Future observations to find M32-analog systems as well 
as simulations of spheroidal systems with similar SFHs to 
M32 may shed light on this issue. 

On the other hand, the inferred SFH for F2 shows that 
the stellar populations of the inner regions of the disk and 
spheroidal components of M3 1 are older and more metal-poor 
than M32. Its mass-weighted mean age and metallicity are 
(Age) =9.15 ±1.2 Gyr and ([M/H]) = -0.10 ± O.lOdex, re- 
spectively. F2 has two main components: 65% ±9% of the 
mass composed by a 8-14 Gyr old more metal -poor popula- 
tion and 30% ±7.5% of the mass in more metal-rich stars of 
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5-8 Gyr old. There is a small contribution from stars younger 
to 5 Gyr to the total mass. The inner disk and spheroidal stel- 
lar populations seem to be indistinguishable from those of the 
outer disk and spheroid. Assuming that M3 1 's disk at F2 1 
disk scale length) has a mean age between 5 and 9 Gyr, our 
results are in agreement with inside-out disk formation mod- 
els. But of course, we need more observations and statistics 
to confirm or rule out this suggestion. 

Lastly, while this paper accounts for the SFH history of the 
bulk of M32's mass, it does not offer strong constraints on 
small "tracer" populations that may testify to the present level 
of any very recent star formation in M32, as well as fossil 
remnants that may date to times well before the 10 Gyr mark, 
after which the bulk of M32 stars were formed. Digging 
down to fainter stars in the M32 CMD to detect any very old ( 
> 10 Gyr) MSTO is not possible with HST, or any instrument 
presently under development, given that the M31 background 
prevents observation at lower surface brightness levels where 
the HST angular resolution would be more effective. Instead, 
we believe the best hope of detecting any ancient, metal-poor, 
remnant from the very first stages in the M32 progenitor 
would come from unambiguously detecting RR Lyrae stars 
at higher surface brightness levels than were observed in Fl, 
where the M3 1 contamination is relatively much weaker Our 
own image simulations show that RR Lyraes can be detected 
with HST at substantially smaller radii than the Fl location. 
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Likewise, better constraints on a < 2 Gyr population native 
to M32 may be provided by bright, if rare, tracers that can 
still be recognized at much higher M32 surface brightnesses 
than were observed in the present work. As with RR Lyraes, 
bright young blue MS stars or AGB stars should be detectable 
throughout the body of M32. 
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^ Color and magnitude shifts of the observed CMD at 
which the „„■„ value is reached. 
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APPENDIX 

EFFECT OF BINARIES 

The resuhs presented in this work were obtained assuming a 35% binary fraction in the synthetic CMD. To investigate how 
much this assumption might affect our solution, we have repeated the entire process of deriving the best mean SFH of Fl and 
F2 assuming not only 35% but also 0%, 70% and 100% binary fractions in the synthetic CMD. The mass ratios between the 
components of the binaries were set to be uniformly distributed between 0.5 and 1 . 

Table|4]shows the values of the min reached for Fl and F2 as a function of the assumed binary fraction, and using the BaSTI 
and Padova/Girardi stellar libraries. We can see that for Fl the goodness of fit does not significantly improve when varying the 
binary fraction if we use the stellar library BaSTI. However, Girardi/Padova models finds the best fit to the observed CMD in Fl 
when the fraction of binaries is 35%. We therefore choose this fraction as our baseline model. For F2, the xtmin ^ function 
of binary fraction is nearly constant, regardless the stellar library used. Note that BaSTI library always recovers a better fit, i.e. 
lower xl mill than Girardi/Padova ones for both Fl and F2 observed CMDs. The position in the ((5(color), (5(magnitude)) grid at 
which xi mill is reached for Fl is nearly insensitive to changes in the model binary fraction. This is not the case for F2, which 
reflects the fact that its CMD is deeper than that of Fl. 



Figure 18 shows the comparison of the derived SFHs. The SFR as a function of time for Fl (left panel) and F2 (right panel) 
indicates mat the calculated solution does not change significantly but becomes older as the number of binaries increases in the 
model CMD. This is expected: the larger the number of binaries in a system, the more luminous the effective (that is, observed) 
MS and flie brighter and redder the effective MSTO of its CMD. 
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Fig. 1 8. — Comparison of the SFRs as a function of age for different assumed binary fractions in tlie syntlietic CMD. The left panel shows the results for Fl 
and the right panel shows the results for F2. The solution becomes older as we increase the fraction of binaries in the model. This can be clearly seen in the first 
two bins of the SFR in Fl (left panel), which represent ages of ~ 10 and ~ 4 Gyr, respectively. This reflects the fact that as we increase the number of binaries 
of the model CMD, its effective MS becomes more luminous and its effective MSTO becomes brighter and redder. 



